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Abstract 

The dynamics of driven spin boson model is studied analytically by means of the perturbation 
approach based on a unitary transformation. We gave the analytical expression for the population 
difference and coherence of the two level system. The results show that in the weak driven 
case, the population difference present damped coherent oscillation (single or double frequency) 
and the frequencies depend on the initial state. The coherence exhibit damped oscillation with 
Rabi frequency. When driven field is strong enough, the population difference exhibit undamped 
large-amplitude coherent oscillation. The results easily return to the two extreme cases without 
dissipation or without periodic driven. 
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1. Introduction 

Discrete electronic states of the qubits can be considered as localized two levels sys- 
tems (TLS), which has been related to the potential use as building blocks of prospective 
quantum logic gates |lL2|. Solid-state qubits, including superconducting , IfJ and 



semiconductor! 



10| ones, couple to environmental degrees of freedom that potentially 
lead to dephasing. The development of techniques to protect coherence of carrier tunneling 
in the qubits are key points for a successful implementation of quantum information pro- 
cessing in these systems and the development of qubit based quantum information systems. 

Since the development of local spectroscopy techniques, the atomiclike optical properties 
of semiconductor quantum dots (charge qubit) or superconductive ring (flux qubit) have 
been intensively studied. The group of T. H. Stievater have reported the first observation of 
Rabi oscillations from excitons confined to single GaAsZAlo.3Gao.7As QD L. Besombes 
et al. optically control the charge state of a single QD and coherently manipulate the 
confined wave function exploiting quantum interference and Rabi oscillation phenomena by 
microspectroscopy in individual InGaAs/GaAs QDs[l3]. Using a pulse technique, Pashkin 
et al. coherently mix quantum states and observe quantum coherent oscillation of coupled 
qubits in the vicinity of the co- resonance [l3j] , the spectrum of which reflects interaction 
between the qubits. Zrenner et al. actualize Rabi oscillation by playing an InGaAs quantum 
dot in a photodiode and demonstrate that coherent optical excitations in the quantum 
dot two level system can be converted into deterministic photocurrents[14]. While many- 
body effects of the environment fundamentally change many aspects of Rabi oscillations, 
particularly the lack of saturation or decoherence. Thus, we think over two questions: first, 
how does the interplay or competition between optical driving and dissipation in the TLS 
influence the decoherence? second, can one sustain large-amplitude coherent oscillation in 
driven spin boson model? A good knowledge of the decoherence holds most prominently for 
many applications such as optoelectronic devices in quantum information processing where 
the operation completely relies on the presence of coherence jis|. 

The dynamics of driven spin-boson model 16] also attract theoretical physics attention 
for its widespread applications to various biological, chemical, and physical systems, e.g., ac- 
driven superconducting quantum interference devices, laser-induced isomerization of bistable 
molecules, laser-induced localization of electrons in semiconductor double-well quantum 
structures, or paraelectric resonance. The various communities typically rely on different 
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methods of description. The most direct approaches is portrayal of the time evolution of the 
corresponding reduced density matrix p(t), which starts from the generalized master equa- 
tion (GME) of the TLS. Two popular approximation are either based on the system-bath 
coupling expansion obtained by use of a projector operator method (commonly known as the 



Bloch-Redfield formalism 



17l|). or on the expansion in the coupling expansion in the coupling 



matrix element A (such as a tunnel splitting) by use of (real-time) path-integral methods, 
such as the application of the so termed noninteracting blip approximation (NIBa[i^]). Lud- 
wig Hartmann et al. have enlighten the advantages and disadvantages of two approaches |l9|. 
A special case where the field frequency is comparable to the TLS frequency (resonance or 
near-resonance) is prominent important in experiment, but is more difficult to handle ana- 
lytically. If the system-bath coupling is weak, the traditional optical Bloch equations will 



produce a meaningful result 



20 



2 if ] . However, their model with decoherence rate as a con- 



stant parameter is too simple to make detailed quantitative predictions here 22J. In general, 
to obtain a solution for time-dependent spin-boson problems even numerically is nontrivial 



task|2J 



24|. 



In this paper we study the quantum dynamics of driven spin-boson model, where the 
driven field is near resonance with the TLS. Analytical explicit expressions for popula- 
tion difference {cr z (t)) and coherence {c x (t)) are presented through perturbation treatment 
based on a unitary transformation. The result shows that with increasing amplitude of 
the driven field, the population difference transferred from damped quantum oscillation 
(single-frequency or double-frequency) to large amplitude undamped coherent oscillation. 
Therefore, one can efficiently control quantum coherent dynamic by optical pulse, induce 
and maintain large-amplitude coherent oscillations. The critical condition from damped to 
undamped large-amplitude coherent oscillations is given. The coherence decays to ground 
state with the Rabi frequency in the case of weak driven. We also investigated the depen- 
dence of the population difference (cr z (t)) and coherence (cr x (t)) on the initial state. 

The paper is organized as follows: In sec 2 we introduce the model Hamiltonian for 
driven spin-boson model and solve it in terms of a perturbation treatment based on unitary 
transform. We analyze the population difference and the coherence in different initialization 
and give discussion in sec 3. Finally, the conclusion is given in sec 4. 

2. The model and theory 

We study the driven spin-boson dynamics, where the two level system is linearly coupled 
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to a continuum of harmonic oscillators and is driven by classical microwave field. The system 
under consideration can be modeled by the HamiltonianQ, 0]- 

H{t) = H S + H d (t) + H b + H l (1) 

with 

H s = -Aa x /2 (2) 

H d (t) = e(t)a z (3) 

H b = ^cj k b+b k (4) 

k 

H l = \Y^9k{bt + b k )a z (5) 

k 

where H s is the Hamiltonian of the TLS, H^t) of the external driven field, H b of the bath 
and Hi of bath and system interaction that is responsible for decoherence. Throughout this 
paper we set h = 1. Here cr, are pauli spin matrices, A describe the coupling between the 
two state and e(t) is the external time dependent modulated field. b^(bk) and us k are the 
creation (annihilation) operator and energy of the phonons with the wave vector k. g k is the 
electron-phonon coupling strength. The ohmic bath is characterized by its spectral density: 

J(lo) = — u k ) = 2auj9{uj c — uj) (6) 

k 

where a is the dimensionless coupling constant and 9(x) is the usual step function. 

Firstly we diagonalize the Hamiltonian independent of time H s + H b + H^. Here we apply 



a canonical transformation, H' = exp(s)H exp(— s) with the generator 25]: 



S = £t7&(#-&*)* (7) 
, ^k 

k 

Thus we get the Hamiltonian H' and decompose it into three parts: 

H' =H' + H' 1 + H' 2 (8) 

where 

H'o = ~\v^ x + E u ^ b " - E £j-^(2 - (9) 

fe k k 
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H 'i = 4 E^ 1 ~ $M + b ^ - &v E -&(K - &*) (io) 



2^ jn " - 2 y ^ u k 

k k 



H' 2 = -lAa x (cosh(£-ZUbt-h)-v) 

^ , u k 1 —f u k 

k k 



with 



i 



^exp(-^^) (12) 



(13) 

L} k + T]A 

Obviously, H' can be solved exactly. We denote the ground state of H as \g) = |si)|{0jt}), 
and the lowest excited states as |s2)|{0fc}), l s i)|{lfc}) where \s\) and \s2) are eigenstates of 
Cx (&x\ s i) — &x\s2) — — I S2) ) , means that there are n k phonons for mode k. The 

last term of H' is constant energy and has no effect to the dynamics behavior. Thus, we 
can diagonalize the lowest excited states H' as: 

H = — r)A\g)(g\ + \^E\E)(E\ + terms with higher excited states (14) 

E 



The digitalization is through the following transformations 25 1 



where 



and 



\s 2 )\{0 k }) = J2*(E)\E) (15) 

E 

l«i>|{u}> = £w(^ (i6) 

E 

\E) = x(E)\s 2 }\{0 k }} + "£y k (E)\ Sl )\{l k }) (17) 



<*)=^Y. {E+ ZL^ )l (18) 



E + ~r]A - uj k 



with Vk = r]Agk£,k/uk. E are the diagonalized excitation energy and they are solutions of 
the equation 

E-l V A-Y A = (20) 

A series of experiments have successfully realized coherent control of the qubit by applying 
resonant microwave excitations [3]. The qubit state evolves driven by a time- dependent term 
e mw cos(2n Ft) a z in the Hamiltonian, where F is the microwave frequency and e mw is the 
microwave amplitude. So we focus our attention on a monochromatic field of the form 
e(t) = e cos(Qt)cr z and choose the control field excitation to be resonant with the splitting 
of the tunneling TLS such that it also does not induce transition to higher-lying excitation. 
So Hamiltonian can approximately be describes as: 

H' = ~ V A\g)(g\ + Y,E\E)(E\-sco S (Qt)a z (21) 

E 

After expand with eigenstate \g) and \E) : 
H' = ~ V A\g)(g\ + Y,E\E)(E\ - ~ ^*(E)\g)(E\ exp(iQt) + x(E)\E)(g\ exp(-int)) 

E E 

(22) 

In deriving Eq.22 we have ignored the nonenergy conserving term or counter rotation 
term in the rotation wave approximation 2^| . This is generally a very good approximation, 
especially in the special case when the two states are at resonance or near resonance with 
the incident field A w Q. In the interaction picture, 

H' 0I = -~r)A\g){g\ + Y,E\E){E\ (23) 

E 

Vi(t) = - £ -Y,(x*(E)\g)(E\eMK-lv&-E + m 

E 

+x(E)\E)(g\ exp(-i(~ V A -E + Q)t)) (24) 

The wave function can be written in the form \^(t)) = C\{t)\g) + CE(t)\E), where 
Ci(t) and C^(t) are the probability of find the particle in state \g) and \E) at time t, 
respectively. The corresponding Heisenberg equation is 

^|*(f)> = -iVKt)|*(f)> (25) 



Assumed that the initial state of the model is exp(S')|\I/(0)) = Ci(0)|si)|{0fe}) + 
^2(0) 1^2) I {Ofc}) = Ci(0)\g) + 6*2(0) ^2 E x(E)\E). Using Laplace transformation, we obtain 



Co(P) 



gi(Q) + g 2 (Q)|EE^S 

km- 1 



(26) 



E E-iP 



and 



(c 2 (o) * (p + i(-kA + n)) + iCMD * ~ l E 



In the complex function theory, the sum of E can be simplified: 
x(P) | 2 1 f dE 



(27) 



E-iP 



277/7 (^-h A -E fclT T^)(P-,P) 
1 



(28) 
(29) 



iP+yrjA-uj k 



After Changing zP + |ryA tow + i0 + |27|, C (P) and x*(E)C E (P) can be rewritten to: 

^ - 77A - R(co) + i 7 (w))C 1 (0) - iC 2 (0)§ 



and 



w - fi) (u; - 77 A - R(co) + ij{ui)) - 
i(co-Q)C 2 (0) -iCi(O)! 



(a; - - r?A - P(^) + i 7 (u)) - \ 
where R(co) and 7(0;) denote the real and imaginary parts of J2k ' /( u ~ u k) 



(30) 



(31) 



R{co) = -2a 



co + 77A 



co r 



CO 



CO r 



r/A co + 77 A 



In 



\co\ (co c + 77A) 
^A(w c - a;) 



(32) 



and 



'-/(co) = 2a7ro>- 



( v Af 



: {6{lo) + 6{lo c -co)-1) 



(w + ijA)*^-' ' ~' *' (33) 

Comparing with Bloch equation and Markovian approximation, decoherence rates j(oo) be- 
come frequency dependent themselves. R(co) is determined by the whole of the spectrum 
of frequencies in the spectral density J(oo) of the ohmic bath. According to the correlation 
of real and imaginary parts in complex function theory, 7(0;) also involve the whole of the 
spectrum of frequencies in the spectral density J{co) of the ohmic bath. That is more general 
and physical. 
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Then we inverse Laplace transformation to time parameter space: 

+00 

= exp(^At/2) r exp r iut + + )duJ (34) 

2tt J 



and 



+00 

y^x*(u)C E (t) = ex P(^ At / 2 ) f x*{u)C E {uj)exp{-tujt + + )du. (35) 



— OO 

The expectation value <jj can be expressed as: 

(ailt)) = (^(0)1^(^)1^(0)) = (tf(0)| ex P (-S)ex P (-iH't)a i exp(iH't)exp(S) |tf(0)> (36) 

After straightforward calculation, the electron population difference (a z (t)) and the coher- 
ence {(J x (t)) are gained as follows: 

(a z (t)) = 2Re(C*{t)Y,^*(uj)C E {t) exp(-zOt)) (37) 

E 

and 

(a 1 (i))==^(l-2^ I *(«^xl^(i)). (38) 

E E' 

Rather simple expression for the population and coherence are obtained analytically. 
3. The result and discussion 

Here we only discuss the near resonant case Q = A. uj c is taken as the energy unit. 
Without special indication, the coupling constant of the environment and TLS is taken as 
a = 0.01. 

In the especial case without coupling to the environment, we rotate a z axis around a y 

axis to a x axis, correspondingly a x axis to a z axis, the Hamiltonian turn into the famous 

quantum optics model describing the interaction of a single-mode radiation with a two-level 
r— 1 

atom|26|. As we have known, the dynamics of (<J z (t)) depend on the initial condition. If the 
qubit is at (a x (t = 0)) = 1 eigenvalue, the driven field is turned on, then, using the rotating 
wave approximation, i.e. in the near resonance condition, the population difference (a z (t)) 
will exhibit quantum beats that result from the interference of fast oscillation with field 
frequency Q and slow oscillations with Rabi frequency e. The coherence (a x (t)) oscillates 
with the Rabi frequency. But if the system is initially in the state (a z (t = 0)) = 1, without 
the environment, (<J z (t)) will present cos fit oscillation and (a x (t)) = 0. We also knew that 



in the other especial case lack of the driven field, (<J z (t)) damply oscillate [25(. The question 



is how will this dissipation environment and driven field collectively influence the population 
difference and coherence of TLS. 

In Fig. 1, we plot coherent dynamics of the population difference (a z (t)) in the initial 
state (<J z (t = 0)) = 1, where a driving field is of the same frequency Q — A, but of different 
amplitude. The curves shown in panels (a), (b) and (c) correspond to three external field 
amplitude e = 0, 0.01, 0.1. Fig. 1(d) present population difference (cr z (t)) with e = 0.4 
and a = 0.3. The Fourier spectrums of the population difference are presented in the right 
four corresponding left panels of Fig. 1. As seen from Fig. 1(a), in the absence of driven, 
population difference exhibits damped oscillation. The frequency uiq and damping rate 7 of 
the oscillations is well agree with prevenient result 25| . With weak driving field strength, 
(<j z (t)) decays with beat pattern in Fig. 1(b) of two frequency almost Q — e and Q. Further 
increasing driving field strength over the regime of half-width 7 of the driving frequency Q, 
(a z (t)) damped oscillate with a master frequency Q in Fig. 1(c). Above three case shows that 
in the weak driven case, the system is always damped and this genuine quantum coherence 
oscillation is weakened by friction from the dissipation environment. For sufficiently strong 
amplitude e > 20, population difference damped a little at initial short time, then present 
undamped large-amplitude coherent oscillations with the frequency of the incident field fl 
That is to say, the driven field overcome dissipation and domain the dramatics of the system. 
It is clearly seen that quantum beat patter of the population difference (a z (t)) oscillate only 
appears in the small regime of driving strength, which is depend on the frequency Q — e 
and/or Q + e whether moves out of the spectral width 7 of the incident field Q or not. 

In Fig. 2, the expectation value (a x (t)) plotted as a function of u c t for TLS driven by 
a monochromatic field with the same parameters as in Fig. 1. As seen, in the absence of 
driving, the coherence (a x (t)) reaches asymptotically a value 1, implying that mainly the 
ground state is occupied at long times 28[]. For weak driving amplitude e = 0.01 and 0.1, it is 
visible that (<r x (t)) oscillated tends to 1 and the oscillation frequency is the Rabi frequency, 
which is directly proportional to driven amplitude. While for large driven field strength 
e = 0.4, the coherence (a x (t)) oscillates around zero with a small amplitude. Driven field 
make the transfer probability to left state and right state equal, that is to say the electron 
tunnel between the two states left dot state \L) and right state \R) , so the undamped 
oscillation occurs as seen in Fig. 1(d). The situation of strong driving field can be correspond 
to the case of absence of damping, where the coherence (cr x (t)) is exactly equal to zero. 



S 



Since any superposition \^Sf) = a\sl) + (3\s2) can be prepared in experiment, through 
manipulation of the quantum state is performed by applying microwave pulses fi(t) to the 
gate 29J. Fig. 3 gives the population difference (<J z (t)) for the same parameters as in Fig. 
1, but assume that at time t = the particle is held at superposition (a x (t = 0)) = 1, with 
the bath being in thermal equilibrium. Without driven field, (<J z (t)) = 0. For weak driven, 
despite the dissipation we were able to induce Rabi oscillations by applying microwave 
pulses at the near-resonance frequency of the TLS. The oscillation decay nonexponentially 
and display a clear beating. From the Fourier spectrum, we find two oscillation frequency, 
one is the frequency of the incident field Q — e and the other is the Rabi frequency Q + e, 
that is to say, (<J z (t)) can be fitted by sinet * sin fit. This corresponds to the common 
quantum beat in the quantum optics. The frequencies of population difference oscillation 
are dependent on the initial state. In Fig.3(d), a well-behaved coherent oscillations can 
be observed after adequate long time when the driven field strength is adequately strong 
e = 0.4. The frequency of the oscillations firstly is Q + e, then Q. 

In the same parameters and the initialization with Fig. 3, coherence (a x (t)) are plotted 
in Fig. 4. In the presence of tunneling with a dissipation bath and weak driven, due to the 
dominance of transitions to the continuum modes, we see that the coherence (a x (t)) presents 
damped oscillation and finally stabilizes the ground state as seem from the Fig. 4(a) and 
4(b). The frequency of the oscillations is Rabi frequency. The damped can be look as a 
leakage of energy 30]. The leakage can be suppressed by enhancing of the external field. For 
strong driving field, the long time limit is shift to nonzero by external field. 

The population difference (<J z (t)) will exhibit damping quantum beat, which have be 
observed in superconducting qubit experiment [3J]. Comparing Fig. 1 and Fig. 3, these curves 
show clearly that the initialization of TLS determine the features of the quantum beat. The 
condition to present beat in the initial state of Fig. 1 is more rigour than that of Fig. 3. The 
quantum beat only emerge in a small regime of the driving field amplitude, that is to say the 
driving field amplitude must be in the range of the damping rate 7 of the driven frequency 
Q. While in the initialization of Fig. 3, the quantum beat can always exist, except the 
driven amplitude enough strong to undamped oscillation. The frequency of the beat pattern 
Q and Q — e or /and Q + e in Fig. 1, While in Fig. 3, the frequency of the beat pattern Q — e 
and Q + e. It is helpful to experimenter for the empirical function choose. The difference 
of frequency in Fig.l and Fig. 3 is understood qualitatively as follows. In the initialization 
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(a z (t = 0)) = 1, Hf(t) = e{t)a z is the maximum of Hf(t) but Hi = \ s ^ Jk gk(p~l + bk)cr z 
is zero due to |{?i.fc}) = |{0fc}), so the driven field firstly domain the oscillation of TLS, 
then environment and driven field interplay. At the end the frequency of the oscillation 
give priority to Q, only when Q — e or /and Q + e in the range of decoherence rate 7 of the 
driving frequency Q, quantum beat emerges. While in the initialization (a x (t = 0)) = 1, at 
the very start, Hf(t) = e{t)a z and fZj = ~ ^2 k gk{b~l + bk)u z are zero together, driving and 
bath interact so the frequency of the oscillation are Q — e and Q + e. 

In what follows, we give the transform condition from damped to undamped oscilla- 
tion. From derivation, we can see that much of the behavior of the population difference 
(a z (t)) and coherence (a x (t)) can be learned from a study of the singularities of Cq(uj) and 
Y2e x*(oo)Ce(^)- According the complex-function theory, we analyze the real part of the 
energy denominators, ReD(cc;)= (u — U)(lu — rjA — R(u})) — that is symmetry in e. Fig. 
5 shows ReD(c<j) with different driving field amplitude e = (solid line), 0.1 (dashed), 0.3 
(dotted), 0.5 (dash-dotted). When e = 0, the singularities of Cq(u) and x*{uj)Ce{u) are 
uj = and u = ljq, which are correspond to zero drive and Lamb Shift due to the dissipation. 
co>o is determined by the function u — rjA — R(u) = 0[25J]. Owing to the driving field, two 
singularities are shift. Define the two singularities U\ and uj 2 and assume u)\ < u>2- We find 
that two positive singularities are in two sides of ujq in the case of small external driving 
field amplitude e. With the increasing of e, one pole U\ becomes larger and the other u 2 
becomes smaller, at the end the small singularity come to the negative axis. Furthermore, 
the imaginary parts of the denominator — Q) can be understood as a damping rate. 

Commonly, they depend on the frequency. Because of the interactions with the continuous 
spectral dissipation bath and driving field, the energy of the localized state is altered to new 
value tO\ and 102- The nature of the solution depends critically on whether the energy u\ and 
u>2 are within the band of state uo k . Therefore, the continuous band of state Uk is confined 
to the range < uo k < u> c , so the character of solution determined whether io\ and U2 are 
also within this range. If the values are within the range < U\,U)2 < uj c , the solution has 
an important property: dissipation domain the behavior in the system, after adequate time, 
the coherence oscillation is damped out. The reason is that generally imaginary component 
of denominator is not zero throughout the continuum band of photon, the decay rate 71(0;) 
and 7i(co>) corresponding to uj\ and u 2 are finite. That is to say, the damping rate 7(0*) is 
nonzero if only in the section of < uj < uj c , which is the key point of the result. Of course, 
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if one of the two energy u\ or uj 2 is outside of the band of continuum state, i.e. u\ is negative 
or c^2 > ^c, then the corresponding damping rate of the correspond energy vanishes. (<r z (t)) 
can't decay to zero and always exhibits large amplitude coherent oscillation. The condition 
damped oscillation or not depend on two roots of ReD(u;) in the dissipation regime or not. 
Approximatively estimate when e > 2Q, population difference can maintain undamped long 
coherence oscillation. We illuminate the possibility to induce and maintain large amplitude 
coherent oscillation by applying a resonant control field and present the condition of the 
long time coherence. It is favorable for quantum computation. The result is similar with 
Ref. 20. 

4. Conclusion 

In this paper, we have investigate the quantum optical control dynamics of the driven 
spin boson model by a perturbation treatment based on a unitary transformation. The 
population difference and coherence are obtained explicitly. Our approach is not restricted 
by the form of spectral density and the initialization. Additionally, a simple expression 
for the real part of the energy denominators, (cu — Q)(u — r/A — R(u)) — allows us to 
analyze the transform condition from damped quantum oscillation to undamped oscillation. 
We find that, for weak field, the population difference (a z (t)) presents damping quantum 
beat, while for strong field e > 2Q, (<7 z (t)) preserve large amplitude coherent oscillation 
with the frequency f2. The coherence (<J x (t)) is also studied. Finally, we hope that this 
work will stimulate more experimental and theoretical works in quantum information and 
computation for quantum optical control. 
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FIGURES 



Fig.l. The population difference (a z (t)) versus u c t in the initialization (a z (t = 0)) = f 
with various values of the driving field amplitude, e — 0, 0.01, 0.1, 0.4, correspond 
to left panels (a) to (d). Also shown with the right panels are the Fourier spectrum 
with parameters same as the left panels, the coupling constant of the environment 
and TLS a = 0.01 except (d) a = 0.3. The driving field frequency near-resonance with 
TLS n = A. 

Fig. 2. The coherence (<J x (t)) versus uo c t in the initialization (<r z (t = 0)) = 1 with various 
values of the driving field amplitude, e = (solid line), 0.01 (dashed line), 0.1 (dotted 
line), 0.4 (dot-dashed line). 

Fig. 3. The population difference versus uj c t in the initialization (cr x (t = 0)) = 1 with various 
values of the driving field amplitude, e — 0, 0.01, 0.1, 0.4, correspond to left panels (a) 
to (d). Also shown with the right panels are the Fourier spectrum with parameters 
same as the left panels, the coupling constant of the environment and TLS a = 0.01. 
The driving field frequency near-resonance with TLS Q = A. 

Fig. 4 The coherence {<J x (t)) versus uj c t in the initialization (a x (t = 0)) = 1 with various 
values of the driving field amplitude, (a) e — (solid line), 0.01 (dashed line), 0.1 
(dotted line), (b) 0.3, (c) 0.4. 

Fig. 5 ReD(c<j) vs uo for different driving field amplitude e = 0(solid line), 0.1 (dashed line), 
0.3 (dotted line), 0.5 (dot-dashed line). 
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